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ABSTRACT 

We report the discovery of a new candidate ultraluminous X-ray source (ULX) in the nearby edge- 
on spiral galaxy NGC 891. The source, which has an absorbed flux of F x ~ 1 x 1CP 12 erg s cm -2 
(corresponding to a L x > 10 40 erg s -1 at 9 Mpc), must have begun its outburst in the past 5 years 
as it is not detected in prior X-ray observations between 1986 and 2006. We try empirical fits to the 
XMM-Newton spectrum, finding that the spectrum is fit very well as emission from a hot disk, a cool 
irradiated disk, or blurred reflection from the innermost region of the disk. The simplest physically 
motivated model with an excellent fit is a hot disk around a stellar-mass black hole (a super-Eddington 
outburst), but equally good fits are found for each model. We suggest several follow-up experiments 
that could falsify these models. 

Subject headings: accretion, accretion disks — X-rays: binaries — galaxies: individual (NGC 891) 



> 

OC 



(N 
O 



X 



1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are non-nuclear 
X-ray sources w ith luminosities Lx > 10 39 erg s _1 
(|Fabbianolll989D . These sources are interesting because 
their luminosities exceed the Eddington limit for a IOMq 
black hole, suggesting that they are either "intermedi- 
ate mass" black holes (IMBH s) of M BH ~ 10 2 - 10 4 M Q 
(jColbert fc Mu shotzkv 1999) or stellar- mass black holes 
seen du ring a special time (s uper-Eddington accretion; 
see, e.g. lGladstone et aLll2009l hereaft er GRD09) or at a 
special angle (i.e., non-isotropic Lr: IKing et al.l ,2001). 
In any case, they are important sources for studying 
black hole physics. For a recent review, see lFeng fc Sorial 
(|20ll . 

In this Letter, we report the appearance of a new ULX 
candidate in the nearby isolated, edge-on spiral galaxy 
NGC 891. The gal axy is thought to be a Milky Way ana- 
log in luminosity (Ide Vaucouleurs et al.lll991l ) and color 
(|van der Kruit fc Searldll98lD. and, like the Galaxy , it is 
a barred spiral (e.g. iGarcia-Burillo fc Guelinl 119951) . In 
Section 2, we describe the XMM-Newton detection and 
other observations, in Section 3 we attempt empirical fits 
to the X-ray spectrum, and in Section 4 we discuss our 
results in context. 

Throughout this Lett er, we use a Galactic column of 
N H = 6.5 x 10 20 cur 2 (|Kalberla et all 120051 ) as well as 
the Wilkinson Microwave Anisotropy Probe cosmology 
(H n = 71 km s- 1 Mpc' 1 , fl A = 0.73, and Q m = 0.27; 
ISpergel et all 120071) NGC 891 is abo ut 9-10 Mpc away 
(Tcmple_elall2005; Tu lly et alf2009l ). and we adopt d « 
9 Mpc. 

2. A NEW CANDIDATE ULX 

A new bright source was discovered in a 133 ks XMM- 
Newton exposure of NGC 891 (obsID 0670950101) with 
an aperture-corrected 0.3 — 10 keV EPIC-pn flux of 
Fx = 1.6 x 10~ 12 erg s _1 cm -2 . An inspection of 
prior X-ray observations between 1994 and 2006, in- 



hodgeskl@umich.edu 

1 Department of Astronomy, University of Michigan, Ann Ar- 
bor, MI 48109 



eluding high resolution Chandra and ROSAT observa- 
tions, demonstrates that this source had not yet been de- 
tected at sensitivities down to Fx < 10~ 15 erg s" 1 cm~ 2 
in the 0.3 — 10 keV band. A 2 ks Chandra observa- 
tion (obsID 14376) indicates source coordinates (J2000) 
02:22:33.45±0.03 +42:20:26.8±0. 5 after process ing with 
the sub-pixel EDSER algorithm (|Li et all 12001 . 

Flux monitoring with the Swift XRT (target ID 35869) 
indicates an overall decay with some hint of variability 
(Figure [2]) , but further monitoring is necessary to estab- 
lish a decay. The Swift values include a minor correc- 
tion for contamination by a nearby X-ray binary with a 
separation of about the Swift point-spread function, but 
Poisson noise is the dominant source of uncertainty. 

Follow-up optical observations, including concurrent 
XMM-Newton optical monitor and Swift UVOT expo- 
sures as well as a snapshot with the Mayall 4m tele- 
scope at Kitt Peak National Observatory, detect no opti- 
cal counterpart to a limiting magnitude V ~ 20 mag. 
NGC 891 is regularly monitored for supernova candi- 
dates, with no detections in the past decade and regular 
exposures in the past six months (A. Filippenko, private 
communication). Likewise, a follow-up 5 GHz observa- 
tion with the EVLA at a sensitivity of 15 /LtJy hm _1 (for 
a beamsize 6 ~ 10 arcsec) found no counterpart. Diffuse 
emission from NGC 891 with a flux of F v « 80 /LtJy bm _1 
is detected at this position, meaning any future attempt 
requires higher resolution. 

Spectral fits (Section 3) indicate an absorbing column 
several times higher than the Galactic value, suggesting 
an extragalactic origin. The source is coincident with the 
disk of NGC 891, but it may be a background quasar or 
blazar. Such a source detected in the X-rays ought to be 
associated with a bright optical/UV counterpart with a 
flux dens ity within an orde r of magnitude of the X-rays 
(see, e.g. IShang et al.ll201ll ). and we would expect rapid 
variability in the Swift monitoring for a blazar. The ab- 
sence of these signatures suggests a physical association 
with NGC 891, but this is not yet certain. 

If the source is in NGC 891, the 0.3-10 keV fluxes indi- 
cate (absorbed) luminosities between 5— 15x 10 39 erg s~ x , 
suggesting a ULX or supernova. Since we have seen no 
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Fig. 1. — 0.3 — 10 keV X-ray maps of the field. Top Left: EP1C-MOS1 image with bright contours of the same (overlaid on other images). 
Top Right: Combined Swift XRT image from our monitoring campaign (21 ks), smoothed by three pixels. Bottom Left: 2003 Chandra 
image (120 ks). Bottom Right: 2011 Chandra snapshot (2 ks). Note the additional new source to the north also seen by XMM-Newton 
has a 0.3 - 10 keV F x ~ 2 X 10~ 13 erg s" 1 cm" 2 . 



radio counterpart, a normal Type la or II supernova is 
disfavored. Although we might also expect a radio coun- 
terpart to a ULX, the radio luminosities expected are 
substantially smaller: using the "fundamental plane" ex- 
pression for radio luminosity a s a function of black hole 
mass and X-ray luminosity in [Merloni ct al. (2003j), we 
expect Lradio < 10 34 erg s _1 for a 10 4 M Q black hole. As- 
suming synchrotron emission near 1 GHz with a spectral 
index a = 0.7, the flux density from such a source is at 
most ~50 /iJy, already below the diffuse radio halo in 
the 10 arcsec EVLA beam. The flux would be even less 
for smaller black holes. The radio non-detection is there- 
fore inconsistent with normal supernovae and consistent 
with a black hole. Therefore, the source appears to be a 
new ULX. Hereafter, we refer to it as NGC 891 ULX1. 
The most recent X-ray observation prior to the XMM- 
Newton exposure places an upper bound of 5 years on 
the outburst. 

Archival Hubble Space Telescope ACS HRC and WFC 
images from 2004 (obsID 9414) reveal a possibly nearby 
source within an arcsecond of the X-ray position, but 
not directly associated (Figure [3]). Based on centroids 
of bright sources in the field, the HST astrometry is un- 
certain to within 0.5 arcsec. The source is faint with 
a FWHM diameter of -0.23 arcsec (the PSF is about 
0.1 arcsec at FWHM), and may be a star cluster. In 
the F555W filter, it has a ST magnitude of ~ 24.3 mag 



(-MF555W ~ — 5-7 m ag, corrected for Galactic Ay ~ 
0.22 mag), and in the F814W filter it is greater than 
23.5 mag (Mp8i4w > —6.4 mag, corrected for Galactic 
Ai ~ 0.12 mag). This source is too dim to be seen in the 
optical data described abo ve. If it is a cluster the ULX 
may be physically related (jZezas fc Fabbianol[200l . 

3. SPECTRAL FITS 

We extracted spectra in the 0.3 — 10 keV band from 
the EPIC-pn and MOS cameras aboard XMM-Newton 
using standard Scientific Analysis System recipes. We 
excluded periods of background flaring during the first 
10 ks and last 30 ks, leaving 92 ks of good time intervals. 
Because the field is crowded, we extracted spectra only 
from within a aperture with r — 25 arcsec (corresponding 
to about 80% encircled energy at 1.5 keV), but reported 
fluxes are aperture corrected. The pn count rate in this 
aperture is ~0.45 cts s , and the task epatplot shows a 
negligible pile-up fraction. The MOS spectra are likewise 
not piled up. 

A two-sided Kolmogorov-Smirnov test finds that the 
light curve (Figure [5]) is incompatible with a constant 
flux (P < 10~ 6 ), but variability is weak, with a fractional 
r.m.s. in 100 s bins of only 13% and no evidence for quasi- 
periodic oscillation in the combined, synchronized light 
curve. 

The spectrum is a featureless continuum aside from the 
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Fig. 2. — Left: Combined synchronized EPIC light curve binned 
to 800 s in GTIs. Right: XMM-Newton, Swift, and Chandra fluxes 
measured since 08/2011 with 90% error bars. 



effects of the absorbing column. In this section , we de- 
scrib e empirical model fits using XSPEC vl2 .7.0 (jArnaudJ 
1996). Each spectrum is binned by 20 counts and we use 
the x 2 goodness-of-fit statistic. An acceptable model is 
one that is not rejected at 95% probability by this metric. 
In comparing models that fail the re gularity conditions 
for the F-test ( Protass ov et al.l 120021) . we calibrate the 
F-test in a method similar to the "p a ramet ric bootstrap- 
ping" described in IProtassov et all (|2002D by using the 
best-fit model (the null model) to generate a large num- 
ber (100-1000) of fake spectra with the XSPEC 'fakeit' 
tool for the same exposure time. We fit each simulated 
spectrum with the null model and the test model to ob- 
tain a distribution of Sx 2 , which we compare to the Sx 2 
value for the data. All of our spectral fits incorporate a 
photoelectric absorption component with Njj frozen to 
the Galactic value, which is not listed with our model pa- 
rameters (Table Q]), as well as an absorptio n component 
that is free to vary. We use TBabs with the lWilms et al.l 
(2000) abundances. 

The simplest models that have been applied to 
ULX spectra are absorbed power-laws and ther- 
mal bremsstrahlung. Although formally a poor fit 
(x 2 /d.o.f= 1860.2/1171), the single power-law model 
(TBabs*POWER,law) with T = 2.24i°o2 is instructive 
in that its residuals show th e high energy "curva ture" of- 
ten associated with ULXs (Rob erts et al.l [20061 ) and ex- 
cess emission below 1 keV (Figure 0J. Since Galactic 
black hole binaries do not exhibit such curvature, ULXs 
are not simply scaled-up versions of ordinary stellar-mass 
black holes. A broken power-law (TBabs*bknpower) 
with a break energy near 4 keV is an adequate fit (Ta- 
ble [1]), but it does not completely remove the excess 
below 0.4 keV (Figure H}. Likewise, a bremsstrahlung 




Fig. 3.— 2004-01-17 HST ACS HRC image of NGC 891 in the 
F555W filter (~5400A). Lighter color is positive emission. The 0.5" 
ULX position error circle is overplotted, corresponding to 25 pc 
across. 



model (TBabs*bremss) is a good fit with kT 3.4 keV 
but retains this excess. These fits can be improved by 
adding a cool thermal (apec) model with fixed kT = 
0.1 keV; Sx 2 — 35 in the broken power-law model for 
one fewer degree of freedom indicates an improvement at 
over 95% significance. In fact, the excess below 0.4 keV is 
present in all our fits (Figure@|, but the pn and MPS dis- 
agree at these energies (see examples in lStuhlinger et all 
2006) so we cannot conclude that the thermal component 
is physical and do not include it in our reported fits. We 
also exclude MOS data below 0.4 keV. However, includ- 
ing the thermal component does not cause large shifts in 
model parameters in Table [TJ 

We attempted some of the disk models recently tried 
on high qual it y XM M-Newton spectra in GRD09 and 
IWalton et all (poll whose source NGC 4715 ULX1 
has a qualitatively similar spectrum). These in- 
clude a multi-colored disk (MCD) blackbody model 
(TBabs*diskpbb), a MCD model with Comptoniza- 
tion (TBabs* (diskpbb+COMPTt)), a disk irradiated by 
Comptonized photons (TBabs*DISKIR), and blurred re- 
flection (TBabs*KDBLUR2*REFLIONX). In the Comp- 
tonization models, the seed photons are assumed to orig- 
inate at the inner edge of the disk. 

For the pure MCD model we use a "p-free" model 
in which T{r) oc r~ p , where p is a free parameter (in 
the normal MCD m odel DISKBB, p is frozen at 0.5; 
iMitsuda et al.l I1984D . The best fit has an inner disk 
temperature T ln as 1.6 keV and p « 0.54, indicating a 
hot, "slim" disk. The temperature can be lowered to 
Fin ~ 1 keV by adding a Comptonization component, 
but based on simulated spectra the additional compo- 
nent does not significantly improve the fit. Because the 
spectrum remains disk-dominated in either case, we find 
a degeneracy between a hot, optically thick corona or a 
cool, optically thin one. GRD09 find that a cool, thick 
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Fig. 4. — The 0.3 — 10 keV EPIC-pn spectrum (binned to 25 counts) with the best-fit "hot disk" model overplotted. Below are residuals 
for each model. Note the rightmost residual for the power- law model is at —10. 



corona is preferred, but we find no global minimum and 
poor constraints due to the weakness of the component. 

We obtain a good fit for a cool disk with a domi- 
nant Comptonization component when the disk struc- 
ture is modified by the absorbing Comptonizcd photons 
(diskir). In this model, the inner disk is irradiated by 
Comptonized photons, many of which are "reflected" . 
Some fraction is absorbed and re-emitted with a quasi- 
thermal spectrum, the reby modifying the disk emission 
(|Gierlihski et al.ll2008l ). Significant modification only oc- 
curs when the Comptonized component is much brighter 
than the disk emission, naturally forcing fits with cool, 
relatively dim disks. Thus, this model is expected to be 
important in the low/hard state — a very different situ- 
ation from that implied by the MCD model. In order 
to fit this spectrum, the Compton hump must occur near 
2 keV, implying a cool corona (kT e ~ 1.4 keV) and a 



cooler disk (Ti n ~ 0.4 keV). The intrinsic absorption is 
also substantially lower than in other models. 

It is also possible that much of the emission is Comp- 
tonized photons "reflected" off the inner edge of the disk 
with substantial line "blur ring" due to re l ativis tic effects 
near the black hole. Like iWalton et al.l 1)20111 ). wc find 
a significant amount of blurring is required such that 
the disk must reach to 1.26i?c with an emissivity index 
<7i n > 7, meaning that the emission is dominated by light 
from the inner disk. 

4. DISCUSSION 

Despite the high quality of the X-ray spectrum, several 
very different models produce excellent fits. Thus, addi- 
tional data at other wavebands or a timing study in the 
X-rays are required to discriminate between them. We 
briefly discuss the implications and predictions of each 
model. 
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A Hot Disk — The hot disk model implies a blac k 
hole mass of less than 4OM (e.g. lSoria fe Kuncidl2008h . 
hence suggesting super-Eddington accretion. Similar 
hot -disk ULXs a r e expl ained as stellar-mass black holes 
by IWinter et a l. ( 2006) , and the recent appearance of 
the source makes sense if, as expected, super-Eddington 
accretion is a transient phenomenon. The Chandra non- 
detection places a prior limit of Lx < 10 37 erg s -1 , which 
is consistent with a low/hard state in a stellar- mass black 
hole but quiescence in an IMBH. The excellent fit with a 
small number of parameters is also a point in favor of the 
pure MCD model, and a disk-dominated spectrum is con- 
sistent with an extreme v ersion of the high/soft state seen 
in Galactic black holes (|McClintock fcRemillardl 120061 
GRD09). 

The unabsorbed model luminosity is Lx = 
2.1 ± 0.3 x 10 40 erg s _1 , significantly higher than 
other disk-dominate d ULXs (jSwartz et alj 12004 120031: 
iRoberts et al.l I2002D and modestly higher tha n the 
brightest hot-disk sources in IWinter et all (f2006h . This 
luminosity exceeds the Eddington limit by 2.5 times 
for a IOOMq black hole. Hence, it does not easily fit 
into the "sequence" proposed by GRD09, in which disk- 
dominated sources are sub-Eddington accretors on the 
extreme tail of the stellar-mass black hole mass distribu- 
tion. As the source maintains the same hardness ratio 
seen in the XMM-Newton spectra during Swift monitor- 
ing, there may be substantial scatter in Lx — Tm plots (cf. 
IMiller et al.ll200l . If the source is thermally dominated 
(the simplest model with an excellent fit ), it would join 
only a few other ULXs seen in this state ([Feng fc Kaaretl 
[2010: 1 Jin et "aT1[20irl ISeFvillat et alJl2011h . 

In the case of super-Eddington accretion we expect 
to see powerful outflows. Since the outburst began at 
most 5 years ago, these outflows would be within ~1 pc 
of the source. We might therefore expect to see ab- 
sorption signatures in a high quality RGS spectrum, al- 
though the source geometry is unknown. Detection of 
such an outflow would place a strong constraint on the 
mechanical luminosity, an important p arameter in mod- 
els of super-Eddingto n accretion (e.g. iBegelmanl 120021 : 
IDotan fc Sha^[20ll . If super-Eddington accretion in 
this system is episodic, broad or multiple absorption fea- 
tures would be expected. 

A Cool Disk — If the hot disk model is similar to 
an extreme version of the high state seen in Galactic 
black holes, a c ool disk may correspo nd roughly to the 
low/hard state (|Gierlihski et al.ll2008[ ). In the irradiated 
disk model, the spectrum is dominated by a cool, Comp- 
tonized component with the subjacent disk component 
dominant at low energies. This model presumes strong 
reflection, but the absence of a strong Fe Ka line is con- 
sistent with the low corona temperature (kT e ps 1.4 keV). 
The emission comes almost exclusively from the very in- 
ner region of the disk, which itself is cool (T- m < 0.4 keV). 
If this disk extends to the innermost stable circular or- 
bit, its temperature would indicate an IMBH (and we 
would expect the luminosity to vary as Lx oc Tj n in fu- 
ture observations). However, the disk may be truncated 
at large radi i outside of a large c orona or outflow (further 
discussed in lFeng fc Sorial 120111 ). in which case the disk 
temperature cannot be used to infer a mass. 
There are a few objections to this model. First, a cool 



TABLE 1 

Spectral Models 



Component 


Parameter 


Units 


Value 


TBabs*bremss 


TBABS 
BREMSS 
X 2 (d.o.f.) 


kT 


10 21 cm -2 
keV 


1.3 ±0.1 

3.4 ±0.1 
1189.3 (1171) 




TBabs* 


BKNPOWER 




TBABS 
BKNPOWER 

X 2 (d.o.f.) 


N h 

r 2 

BreakE 


10 21 cm -2 
keV 


1.9 ±0.1 
1.91 ±0.04 
3.2 ±0.2 

o 7 +0.2 

1217.0 (1169) 


TBabs*diskpbb (Hot Disk) 


TBABS 
DISKPBB 

X 2 (d.o.f.) 


N H 
T in 
P 


10 21 cm" 2 
keV 


1.5 ±0.1 

1 62+ 05 

1 - DZ -0.06 

0.54 ±0.01 
1169.1 (1170) 


TBabs*(diskpbb4 


-COMptt) (Hot D 


isk) 


TBABS 
DISKPBB 

COMPTT 

or 

X 2 (d.o.f.) 


N H 
T in 
P 

kT e 

Tp 

kT e 

Tp 


10 21 cm- 2 
keV 

keV 
keV 


1.5 ±0.1 
1.1 ±0.1 
0.54 (f) 
44 a 

0.01 ab 
20 ab 

5.4 a 

1166.6 (1168) 


TBabs*diskir (Cool Disk) 


TBABS 
DISKIR 

X 2 (d.o.f.) 


N H 
T m 

r 

L c /L d 
kT e 

./in 
/out 

r irr 


10 21 cm" 2 
keV 

keV 

fin 
Tin 


o.38t8:°t 

-, Q o+0.05b 
1 - 9c> -0.02 
7 +4 
' -2 

1.4 ±0.3 
0.1 (f) 

04 +0.06b 

g+3 
°-2 
10 5 (f) 

1160.8 (1166) 


TBabs*kdblur2*reflionX (Blurred Reflection) 



TBABS 


N H 


10 21 cm" 2 


2.1 ±0.2 


KDBLUR2 


Rin 


Ra 


L2 6+0.08b 




Rout 


Rg 


400 (f) 




i 


deg 


30 ±11 




9in 


8+1 










<?out 




3.0 (f) 




^brcak 


Rg 


20.0 (f) 


REFLIONX 


r 




1.73 ±0.06 




t 


erg cm s" 1 


> 6000 b 




A Fc 




8 ± 3 


X 2 (d.o.f.) 






1242.34 (1167) 



Note. — All fits incorporate a TBabs component frozen at the 
Galactic Nh — 6.5 x 10 20 cm -2 . Additional absorption is listed here in 
units of 10 21 cm -2 . Errors are quoted at the 90% confidence interval 
based on the XSPEC task steppar. 

a These parameters are not well constrained and we do not quote errors. 
See text. 

These parameters are near the boundary of parameter space. 
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disk component is only a good fit when irradiated, oth- 
erwise a hot disk is required (even the "cooler" hot disk 
has Tj n ~ 1 keV) . If the disk were intrinsically cool, we 
would expect a non-irradiated cool disk to at least be 
competitive with the hot disk. Second, if the black hole 
is in the low/hard state, the prior X-ray non-detections 
suggest it was previously quiescent; even the naked disk 
in this model should have been detected. 

Fortunately, this model can be falsified in a few 
ways. The intrinsic absorption in this fit is small 
compared to the others, so if an optical counterpart 
could be established and its intrinsic color determined, 
a column could be measured. We would also ex- 
pect to see variability like that in Galactic low/hard 
st ates (and potentially quasi-per iodic oscillations as 
in iStrohmaver fc M ushotzkvl 120091 ) , so long-term X-ray 
monitoring can potentially discriminate between a high 
and low state. 

Blurred Reflection — In contrast to the other models, the 
blurred reflection model fits the spectrum assuming it is 
dominated by reflection of Comptonized photons near 
the inner edge of the disk, with relativistic effects smear- 
ing out the strong emission lines. This model predicts a 
Compton hump in the 10-30 keV range, so observations 
above 10 keV could easily distinguish between this model 
and the others. Whereas many ULXs are in galaxies 
with nuclear sources of comparable or greater brightness, 
NGC 891 ULX1 is presently the brightest X-ray source 
in the galaxy. As in the irradiated disk model, blurred 
reflection suggests the low/hard state. The caveats of 
the physical i n terpre tation here are discussed in detail in 
iWalton et ail pOlTI ). but we note here that this model 
requires suppression of the disk emission. 

Like other ULX spectra, the mostly featureless contin- 



uum of NGC 891 ULX1 admits several different scenar- 
ios, suggesting observations in other wavebands or in the 
time domain are necessary to untangle the origin of the 
X-rays. Given that the outburst is recent and preceded 
by at least a 20 year lull, a search for an ionized neb- 
ula like those seen arou nd other nearby, bright ULXs 
(|Pakull fc Mirionl [2002) is of interest. Such a nebula 
would point to recurrent ultraluminous activity, and its 
size and power may allow measurement of the duty cycle. 

5. SUMMARY 

A recent XMM-Newton observation revealed a bright 
new source towards NGC 891. At this point, the most 
natural explanation is that it is a source within the 
galaxy, in which case its flux and the absence of detec- 
tions in other wavelengths makes a ULX the best ex- 
planation. The XMM-Newton spectrum is morpholog- 
ically similar to disk-dominated sources in the GRD09 
"sequence" of ULXs and is indeed fit well by a hot MCD 
disk. However, several very different physical scenarios 
produce good fits to the spectrum, including an irra- 
diated cool disk with strong but cool Comptonization 
and a blurred reflection model. In the hot-disk model, 
NGC 891 ULX1 has a mass less than 4OM and there- 
fore accretion rates over 5 times Eddington. If so, the 
recent ignition suggests a search for outflows would be a 
worthwhile test of super-Eddington models. 

The authors thank the referee for catching errors and 
making comments that significantly improved the pa- 
per, as well as R. Mushotzky for helpful comments and 
A. Filippenko for information regarding SNe monitor- 
ing. EHK gratefully acknowledges support from NASA 
ADAP grant #061951. 

Facilities: XMM, EVLA, Swift (XRT), CXO (ACIS), 
Mayall, HST (ACS) 
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